Mesophases of soft-sphere aggregates 
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Soft spheres interacting via a hard core and range of attractive and repulsive "soft-shoulder" po- 
tentials self-assemble into clusters forming a variety of mesophases. We combine a mean field theory 
developed from a lattice model with a level surface analysis of the periodic structures of soft-sphere 
aggregates to study stable morphologies for all clustering potentials. We develop a systematic ap- 
proach to the thermodynamics of mesophase assembly in the low-temperature, strong-segregation 
and predict a generic sequence of phases including lamella, hexagonal-columnar and body-center 
cubic phases, as well as the associated inverse structures. We discuss the finite-temperature correc- 
tions to strong segregation theory in terms of Sommerfeld-like expansion and how these corrections 
affect the thermodynamic stability of bicontinuous mesophase structures, such as gyroid. Finally, 
we explore the opposite limit of weakly-segregated particles, and predict the generic stability of a 
bicontinuous cluster morphology within the mean-field phase diagram. 
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I. INTRODUCTION 

New classes of soft materials assemble into a bewilder- 
ing array of complex mesophases. These include bicon- 
tinuous structures with multiple symmetries [l|, unusual 
lattices such as the A15 0, well as the more ex- 

otic Frank-Kasper and quasicrystalline packings of dcn- 
drimers and block copolymers [E,i6„i] • This rich behavior 
holds out the promise of designing materials to achieve 
a complex target structure on length scales long enough 
to be used as optical devices or as components of hybrid- 
photovoltaic structures. Nevertheless, the principles con- 
necting structure and interactions at the molecular scale 
to the structure and thermodynamics of long-range as- 
semblies remain poorly understood. 

In this article, wc study the stable morphologies of 
a class of interacting particles with a hard core and an 
isotropic, soft corona. The class is exemplified by a 
corona which costs a fixed energy e when it overlaps with 
the corona of another particle. At sufficient densities, this 
interaction results in stable, finite size clusters arranged 
into a number of structures [1, [HI, [13, [H, [3, 

[lH, [l^. In two dimensions, these clusters arrange into ei- 
ther hexagonally-packed clusters, stripes or a dense phase 
of hexagonally-packed voids d, [1, [l^, [13, [13 ■ I'^ three 
dimensions, numerical calculations of the ground state 
show a number of distinct lamellar, cylindrical and spher- 
ical clusters as well as their inverse form-like structures 
[3, [13, [HI- This behavior is by no means specialized, how- 
ever, as a range of particle interaction types have demon- 
strated to mesophase clustering. Another classic example 
of cluster morphologies occurs for potentials with long- 
range repulsion and short-range attractions [l3, [13, [loj ■ 
Cluster phases have also been demonstrated with mag- 
netic dipoles in a magnetic field in two dimensions both in 
simulations (2]| and, more recently, in experiments (2^ . 

Beyond isotropic particle systems, it is striking that 
predicted cluster morphologies bear a strong resemblance 



to the phase morphologies of block copolymer melts [23| . 
Despite the obvious differences in the microscopic struc- 
ture of polymeric amphiphiles and isotropic particles 
with hard cores, we find that the gross features of the 
long-range order expected in both systems to be in exact 
correspondence both in terms of topology of aggregate 
structures as well as the symmetry of their periodic ar- 
rangement. 

In this paper we explore the full potential of isotropic, 
clustering particle systems to form long-range ordered 
mesophases and investigate the thermodynamic princi- 
ples underlying their assembly. In further analogy to the 
theory of block copolymer thermodynamics, we develop 
a mean-field approach to clustering morphologies which 
is tractable both in the limit of strong-segregation [2^ . 
valid at low temperature and sharp cluster boundaries, 
and weak-segregation [23l, valid near the order-disorder 
transition. In addition to predicting the ground state 
morphologies, our low-tcmpcrature expansion allows to 
probe the explicit relationship between particle entropy, 
interactions and cluster geometry. This expansion high- 
lights a frustration between the particle interaction po- 
tentials and suggests a mechanism whereby cluster mor- 
phologies with high volumetric surface areas are stabi- 
lized at higher temperatures. Interestingly, we predict a 
stable orthorhombic, Fddd bicontinuous phase in a nar- 
row region near the order-disorder transition in coniplcte 
analogy to the prediction [13, [13| and observation [H. [28j 
of this same structure in copolymer melts. 

Previous studies of clustering from isotropic particles 
have largely focussed on properties of a specific type of 
interaction potential. However, in ref. |16| a more gen- 
eral classification of interaction potentials is introduced. 
Here, we take this one step further, demonstrating that 
the gross properties of equilibrium assembly are rather 
insensitive to the precise details of the potential. This 
universality stems largely from universal geometric con- 
siderations both at high temperatures, where this is ex- 
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pected (29[, and at low temperatures. 

In section |lTJ we discuss the general features of soft 
assemblies from the point of view of the Fourier transform 
of the interaction potential. In section Hill we derive the 
phase diagram using a strong-segregation theory, valid 
when the cluster interfaces are sharp and the temperature 
is small. In section ITVl we expand on these results with 
a finite temperature expansion. Finally in section |Vl we 
present a phase diagram valid near the order-disorder 
transition. In section IVll we summarize our results. We 
also discuss a paradigm for designing soft matter clusters. 



II. STATISTICAL MECHANICS OF CLUSTER 
MESOPHASES 

A. Clustering instability and particle interactions 

We begin our analysis by dividing the interaction be- 
tween particles into two parts, the hard-core VHcir) po- 
tential, where 



Vncir) 



oo, 
0, 



r < a 
r > a 



(1) 



and "soft-shoulder" Vss (r) which characterize the longer 
range behavior of particle interactions. On its own, the 
ground states of the step- function corona described in the 
introduction 



Vssir) 



e, r < a 
0, r > a 



(2) 



solve a particularly simple optimization problem. Since 
this potential is purely repulsive, the different morpholo- 
gies that are observed, apparently, minimize the number 
of overlaps within a distance less than a at fixed den- 
sity [l^. It is not clear, however, how other types of 
clustering systems can be explained with such a simple 
geometrical optimization. Moreover, designing clustering 
particles that exploit 'soft packing' requires an additional 
understanding of what principles guide cluster morphol- 
ogy. 

By examining the soft-shoulder potential from the 
point of view of its Fourier transform, 



(3) 



the more general conditions for the clustering behav- 
ior become apparent. In refs. [lB| and [l^, the au- 
thors observed that, at least within a mean- field ap- 
proximation, interaction potentials whose Fourier trans- 
forms have a negative region exhibit an instability lead- 
ing to clustering. Precisely, the instability occurs when 
1 -I- l3V{q)p{l — p/pm) < 0, where p„i is the maximum 
packing density of cores, /3 = l/ksT, p is the average 
density p^ . 

It is, perhaps, not surprising that certain combinations 
of attractive and repulsive interactions can satisfy this 




FIG. 1: (Color online) Two real-space interparticle soft- 
should potentials Vss(t) are shown in (a): eq. ([2]) in blue 
and eq. Q in red with A = 2e, B = e, Aa = cr/3 and 
As = 2(t/3. In (b) the Fourier transforms of each poten- 
tial, normalized by their respective values at g = 0. Both the 
long-range repulsive/short-range attractive and step repulsive 
potential exhibit a minimum for which V{q*) < 0, indicating 
an instability to clustering at sufficiently low temperatures. 



condition. For example, consider the interaction poten- 
tial 

Vssir) - -^cxp (-rVA^) + Bcxp (-rVA|) . (4) 

where A and B are positive and Xa/^b < 1. There- 
fore, at distances r > Xb the soft-shoulder potential is 
repulsive but when r < A/i it is attractive. The Fourier 
transform is 



v^(q) 



^^^.cxp(-AigV4) + ^^'^ 



2\/2 



2^2 



exp I 



X2 



(5) 

which has a single negative peak at a non-zero wave vec- 
tor, provided that AX\ < BX%. More interesting is that 
purely repulsive potentials, like the step potential in eq. 

can exhibit the same negative peak. Figure [T] shows 
that these two very different soft shoulder forms give rise 
to a similar properties of V{c{) in Fourier space, namely 
a minimum at some finite wavevector |q| = q* for which 
l^(q*) < 0. Not even sharp discontinuities are required: 
simulations (2l| and experiments in two-dimensions 
of a dipolar magnetic repulsion softened by an attractive 
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component show clustering behavior even in regimes in 
which the total potential is repulsive. 

From this point of view the fundamental frustration oc- 
curring in clustering particle systems becomes transpar- 
ent. The soft shoulder portion of the interaction poten- 
tial drives the system to order at a particular wavelength, 
A* = 27r/g*, associated with minimum in V{q). At the 
onset of order, this may be accomplished by a weak sinu- 
soidal modulation of the density profile. At lower tem- 
peratures, the soft-shoulder interaction prefers a much 
stronger density modulation, favoring particle packings 
that are at higher densities than can be accommodated 
by the hard-core of the particle. Hence, the hard-core 
portion of the interaction prohibits "over-crowding" of 
particles, forcing aggregates to spread out into clusters 
of various morphologies. 

This observation provides an alternate methodology 
for soft sphere packing problems involving the features of 
the Fourier transform of the potential. Moreover, similar- 
ities in the Fourier transform of two potentials of very dif- 
ferent origins lead one to an effective optimization prob- 
lem expressed entirely in terms of minimizing the number 
of corona overlaps subject to constraints on the particle 
density. 

In further sections of this article, we will focus on po- 
tentials whose Fourier transforms-like those illustrated in 
Fig. ([Ij-are dominated by two modes: the modes at q = 
and |q| = q*. This requires that \Viq)/V{q*)\ < 1, 
for wave numbers at other negative peaks q ^ q*, as we 
find for the potentials discussed in this section Ad- 
ditionally, we assume that V^(q) is maximum and finite 
at q = 0. 



B. Lattice model particle clustering 

To strip away the complexity arising from the hard 
sphere packing of the particles, we follow ref. and 
introduce a lattice gas model to describe the mesophase 
behavior of hard core/soft shoulder particles. In or- 
der to prevent overlaps of particle cores, this model as- 
sumes particles to occupy a close-packed lattice of cores. 
Though this type of model only approximates the trans- 
lational entropy of particle cores, we expect the model 
to become accurate in the low temperature regime where 
clusters are densely packed and when core size is much 
smaller than dimensions of the clusters. A* 3> a. The 
state of each lattice site i is represented by n,; = or 
1. The soft-shoulder pair interaction between sites i and 
j is represented by Vgg = Vssdxj — Xj\). The effective 



Hamiltonian is then written as 



(6) 



where is the chemical potential that controls the den- 
sity of particles in the system. After a delicate Hubbard- 
Stratonovich transformation, the grand partition func- 
tion can be written in terms of a coarse grained field 



0(x) m by Z = J'D(j){K) exp {- A[(j){K)]} , where the ac- 
tion A is 



p„ / (fx In [1 + exp(/?/Lt + i(/i(x))] , (7) 



where, 



d^q 



(27r)3 y(q) 



(8) 



The convergence of the path integral requires that (jj be 
purely real for Fourier modes where V^(q) > and purely 
imaginary for modes where V^(q) < [s^l- The mean- 
field analysis we pursue here is insensitive to this techni- 
cality. 

We pursue our analysis of cq. ([7]) by making two 
important approximations. First, we specialize to the 
mean-field limit, in which (/'(x) is given by the critical 
point of 74[(/)(x)]. Minimizing the action with respect to 
(/)(x) leads to the mean-field equation, for each Fourier 
component. 



(q) = i/3V{q)p{q) , 



(9) 



where pq = V^^ J d?x p(x)exp(— iq • x) and V is the 
total volume. Along with the mean-field condition we 
also have the relation between 0(x) and mean particle 
density for the lattice gas model 



Pn 



g/3Ai+i0(xi) 
1 -)_ g/3M+#(xi) ' 



(10) 



Together eqs. (jU and pO)) form a self-consistent set of 
equations for the mean-field density profile, p(x), and 
self-consistent field, (/)(x). 

The validity of this saddle-point approximation is de- 
termined by strength of interparticle interactions relative 
to the thermal energy scale. Assuming that mean-field 
solutions c>c il3V{(\ = 0) and the mean-field action 
becomes linear in /3V'(0)p^V. The mean-field analysis, 
therefore, requires a steep saddle, or sufficiently strong 
particle interactions l3V{Q)p„i ^ 1. 



C. Level set representation of mean-field solutions 

Clustering behavior is observed when the Fourier 
transformed potential ^(q) takes on negative values for 
some q. This suggests that equilibrium field configura- 
tions will be periodic structures whose fundamental mode 
will have length |q| = q* . Further, from the mean-field 
condition eq. ([5]), we see that if V^(q) is dominated by 
q = Q and q = q* modes, we may consider configura- 
tions of 0(x) that contain only these wavevectors. Note 
that this applies equally to the case of strong-segregation 
where the mean-field density has many non-zero Fourier 
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phases 


G 


m 




ranges of level set tpo 


lamellar 


±x 


2 


2 cos(a;) 


-2 < Vo < 2 


hexagonal 


±x,±\x± 


6 


2 [cos(2:) + 2 cos(f ) cos(^y)] 
L J 


-2 < 7/;o < 6 


sc or P surface 


±x, ±y, ±2 


6 


2 [cos(a;) + cos{y) + cos(2)] 


-6 < i/^o < 6 


bcc or IWP surface 


±x ± y, ±y ± z, ±z ± X 


12 


4 [eos(a;) cos(j/) + cos{y) cos{z) + cos(a;) cos(2)] 


-4 < Vo < 12 


fee 


zbi ± y zb 2 


8 


8 [cos(a::) cos(j/) cos(2)] 


-8<ipo <8 


gyroid 


±i ± y, ±y ± z, ±2 ± X 


12 


4 [sin(y) cos(2) + sin(2) cos(a;) + sin(a;) cos(y)] 


-6 < i/;,) < 6 



TABLE I: Level surface functions for various phases. G is the primary reciprocal vector and m is the number of modes for 
|G| = q* . Notice that the gyroid phase is constructed by a phase shift Sa = n/2. 




FIG. 2: (Color online) Level surface morphologies in terms of level set constants '00 for ImSm space group. I-WP surface with 
^0 = 1-0 (a) starts to pinch off at tpo = 0.0 (b) and evolves to bcc phases ipo = 0.4 (c), tpo = 4.0 (d), and tpo = 10.0 (e). 



modes beyond its fundamental mode. Thus, for a given 
periodic structure, characterized by a set of reciprocal 
lattice vectors, G, we write 0(x) as a superposition of 
plane waves with only q = and q ~ q* modes: 



we may write the mean-field solution form of the density 



\G\=q' 



(11) 



For a given symmetry group, we can take advantage of 
the fact that Fourier transforms along different reciprocal 
lattice directions have the same amplitude to rewrite field 
solutions as 



(x) - 00 = (AgV'(x) 



where 



\G\=q' 



(12) 



(13) 



where Sc, are mode-dependent phases. For a given sym- 
metry, all information about the three-dimensional, equi- 
librium cluster morphology is encoded in iJj{x), the level 
set function. 

Using the mean-field values of ipo and V^g from eq. (O 
and defining, 



P*\V^\ 



and 



(14) 



(15) 



as, 



P(x)/pr, 



ga('i/i(x)-i/'o) 
1 ga(V'(x)-j/'o) ■ 



(16) 



Here, we use the notation Vq = V{0), — V{q*), po = 
/5(0) and = p{q*) ■ Regions of high (low) particle 
density correspond to regions where ip{x) > ipo (-(/'(x) < 
-00); the boundaries of the clusters are defined by the 
surfaces where '0(x) = i/jq. Hence, adjusting the level-set 
parameter ■0o adjusts the size, structure and surface of 
aggregates |31| . 

In table [H we list level set functions which correspond 
to single-mode plane- wave superpositions of the type de- 
fined in eq. (|13p . These include the one-dimensional 
lamellar structure and two-dimensional hexagonal colum- 
nar structure. Three-dimensional structures include cu- 
bic arrays of spheres: simple-cubic (sc); body-centered 
cubic (bcc); and face-centered cubic (fee). Addition- 
ally, these include multiply-connected bicontinuous struc- 
tures, which we refer to by the name of the minimal sur- 
face of the same topology: IWP surface; Schwarz's P 
surface, otherwise known as the "plumbers nightmare" ; 
and the gyroid (G) surface [3l|. In the representation 
of equation I13[ two distinct morphologies can share the 
same level surface function '0(x). For example, the bicon- 
tinuous I-WP surface, having ImSm space group, occurs 
in the parameter range —4 < -00 < 0, whereas the bcc 
structure occurs in < -00 < 12- Similarly, Schwarz's P 
surface shares the same level function as the sc phases: 
sc phases occur in the range 2 < j^ol < 6 whereas the P 
surface occurs when —2 < -00 < 2 (PmSm space group). 
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In Fig. [21 we display a representative example of cluster 
shapes with increasing ipQ for the Imim space group. 

In terms of the mean-field solutions, themselves de- 
scribed by a family of level-set functions, we now have 
the grand potential in the saddle-point approximation, 

^ = f (HKIP^-Vbpg) (17) 

- Pm j d^x In [1 + exp{a{ip{x) - -0o)}] , 

where m is the number of plane- wave modes in iplx). 

III. LOW TEMPERATURE, 
STRONG-SEGREGATION LIMIT 

In the limit of low temperatures the parameter a 3> 1, 
indicating from eq. pB]) that clusters are separated by 
sharply defined boundaries dividing regions of densely 
packed clusters with /o(x) = p,„ from empty regions with 
p{x) = 0. In absolute limit of a — > oo, the temperature 
dependence of all quantities is eliminated. In particular, 
the value of the grand potential for a particular mean- 
field solution has the form, 




FIG. 3: (Color online) Zero temperature phase diagram in 
strong segregation limit. Lamellar, hexagonal-columnar, bcc 
and their inverse phases (from the holes) are found as stable 
structures. 



»(|K|,m) 1 / If. I .2 

-IK 1/5* / ^(^(x)-^o) , (18) 

J in ^ 

where the scaled variables are defined by J7 = ri/Vb, K = 
14 /Vq, and fi = ^/Vq, and the integral is carried out only 
inside the clusters. 

We compute J7 from numerical calculations of po a-nd 
in terms of the level-set parameter for candidate mor- 
phologies. Here, we are assisted by the fact that p(x) 
is nearly p„i in the cluster interior, where ip > ipo, and 
is nearly between clusters, when ip < ipQ. Due to the 
non- linear relationship implied by cq. pop . this step- 
wise density modulation occurs at low-temperature even 
when (f>{'x.) is dominated by a single mode. For lamel- 
lar phases, we can compute the dependence of structural 
quantities on ^pQ analytically. We find po = cos~^{iPq)/tt 
and /5* = sin[cos~^ (?/'o)]/7r. 

We construct the phase diagram as a function of 
(|K|, /i) comparing the computed values of grand poten- 
tial from the candidate periodic phases as well as from the 
uniform state with = 0. Our computed phase diagram 
is shown for the range \Vt\ < 1 in Fig. [31 We find that 
lamellar, hexagonal-columnar, bcc phases, and their in- 
verse phases are stable phases in the strong-segregation 
limit. All transitions between phases of different sym- 
metry are first-order. Note that the diagram is sym- 
metric with respect to p/Vq = 1/2, which arises from 
the "particle-hole" symmetry of the lattice Hamiltonian. 
The relation Eq. (??) returns ji' = 1 — fi for the holes by 
Po = 1 - Po, P'* = "P* and = ~tpo- 
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FIG. 4: (Color online) The grand potential differences from 
the lamellar phase for various candidate structures as a func- 
tion of the chemical potential /i/Vo. 



In Fig. m we also plot the grand potential energy dif- 
ferences from lamellar phase as a function of p for various 
structures along the line of = 0.2. Interestingly, the 
gyroid phase is the second most stable phases in the range 
of 0.4 < /t ^ 0.6 and the most stable among the bicon- 
tinuous phases by far. It is more stable than the lamellar 
phase at very small and large p, though it is pre-empted 
by the hexagonal phase. 

Though our model is unable to distinguish between 
structures that differ only by how the hard cores pack. 
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our phase diagram is in qualitative agreement with sim- 
ulations and a direct minimization of the repulsive- 
shoulder energy that allows lattice packing of cores to 
adjust [ll|. Moreover, our analysis suggests that this 
phase diagram is universal among particle interactions 
whose Fourier transform has a sharp, negative peak. 

IV. FINITE-TEMPERATURE CORRECTIONS 

At finite temperature, but still in the mean-field 
regime, the clustering thermodynamics are modified by 
the corrections arising from the finite width of the in- 
terface between the high and low density regions. This 
intcrfacial width is given by A = l/(a|VV'|) and depends 
on both the large parameter a and the interfacial geom- 
etry. This interfacial region at the boundary of clusters, 
is the primary source of particle fluctuations in the low 
temperature limit, giving rise to entropic contributions 
to the zero-temperature theory described in the previous 
section. 

To compute the corrections to the energy, we employ 
a Sommerfeld-like expansion. It is most convenient to 
compute our corrections from an alternate expression for 
the grand potential 

f3n=^J d^xd^x'p{^)Vss{^ - x')p(x') + 

J d^xS{^)-Pl2po, (19) 

where 

5(x) = pln{p/p,n) + {p„i - p) ln(l - p/pto) (20) 

is local contribution to the entropy of a lattice gas. This 
form is equivalent to the standard density-functional the- 
ory and can be derived from the grand potential using the 
mean- field equations. 

Both inside and outside of clusters where p(x) is equal 
to pm and respectively, the local contribution to the 
entropy is 5(x) = 0. In the low temperature regime 
5(x) is only non-zero at the cluster interfaces where p,„ > 
p{yi) > 0. We expand the grand potential in powers of 
'(/'(x) — ^jjQ and use a Sommerfeld expansion for po and 
/5*. To lowest order in T we have the low temperature 
correction to the grand potential Eq. 

— = padpa - m\Vt\ptSp^ + TSS - jlSpa - paSp, (21) 
where 

TT^ d'^po 

oa^ dtpQ 
6p = Spo - \V*\{Sp^)ipo 



0,004 




FIG. 5: (Color online) Grand potentials for the gyroid phase 
as temperature increases, = 1/10000,1/2000,1/1000 at 
IV, I = 0.2. The calculation from the lowest order corrections 
(solid lines) are compared to the full mean-field calculation 
(dashed lines). Note that the finite-temperature expansion 
becomes less accurate both as T is increased and p/Vo ap- 
proaches and 1 (vanishing cluster width). 



and T = fcsT/Vo. Note that all properties about the 
low-temperature fluctuations is encoded by the level-set 
geometry, by first or second derivatives of po and p^, with 
respect to the level-set parameter ipQ. This also allows us 
to assign geometrical meanings to these terms. For exam- 
ple, the entropy correction SS is proportional to the intcr- 
facial volume, roughly A / dA. When = const., im- 
plying a constant interfacial width, this correction simply 
scales like the interfacial area. Similarly, d^po/dipQ can be 
associated with the curvature of the interface weighted by 
the interfacial width. Since SS represents the entropy of 
enhanced core fluctuations at the interface, its contribu- 
tion to the grand potential is always negative. This sug- 
gests an entropic preference for cluster morphologies with 
a high volumetric surface area, a negative surface energy. 
This tendency for large surface area structures competes 
with the terms derived from soft-shoulder interactions 
already included within strong-segregation theory. 

As an illustration of the thermodynamics of low- 
temperature corrections described above, we consider 
the analysis at p = 1/2 (or cquivalently, po ~ Pm/'^^ 
or -00 = 0); where the lamellar and gyroid phases arc 
shown to be the most and second most stable structures 
within strong segregation theory. For both structures, 
?/'o = is a point of inflection for the density so Spo = 0. 
The contribution to finite temperacture corrections only 
arises from SS and i5/5*. In this case of the fixed den- 
sity, we generally consider (5/3* as to be negative because 
the system becomes more homogenous as T increases 
and the interface between core and cluster broadens. At 
1 14 1 = 0.2, the corrections to the grand potential per 



7 



0.008 r 



bcc(f 2=1 10000) 




iex(f== 11 00(10) 



hex (7-2 1'2000) 
hcx(f!.l/10flO) 



gyroid(f 2= 1/10000) 
gyroi<i(f2=l'2000) 
g>Toid(7-2 1 1 000 1 



0.1 0.2 0.1 0.4 0.5 0.6 0.7 



FIG. 6: (Color online) The comparison of grand potential 
energy among lamellar, hexagonal-columnar, and bcc phases 
from the full mean-field calculation for temperature — 
1/10000, 1/2000, 1/1000 at |V;| = 0.2. 



unit volume (A^) from p*, ~m\V■^,\p■^,Sp^^,, are evaluated as 
SUi, 



4.1123 T2 and dU„ 



5.2592 T^, for lamellar 



and gyroid phases, respectively. Although the interac- 
tion energy of the gyroid phase increases more rapidly 
with temperature than that of the lamellar phase, we find 
that the mixing entropy correction plays a more dom- 
inant role in the overall grand potential energy change 
(per unit volume), giving SSiam = —8.2247 and 
SSgyr = —10.5197 T^, for lamellar and gyroid phases, 
respectively. The surface area of gyroid phases signifi- 
cantly lowers the grand potential and thus become more 
stable at finite temperature relative to the zero temper- 
ature limit. We plot the grand potential changes for the 
gyroid phases as temperature increases (T^ = 1/10000, 
1/5000 and 1/2000) in Fig. (5] Due to the invariance 
of the mixing entropy S under inversion (p p„i — p), 
the contribution of temperature corrections is symmetric 
with respect to p = 1/2. 

This lowest order calculation is expected to break down 
where the interfacial width A becomes larger than the ac- 
tual cluster size of order A* . This can occur either when 
temperature becomes large, or when p* becomes suffi- 
ciently small. This latter condition is always met as the 
density of particles or holes goes to zero at the bound- 
aries of the segregation regime. Therefore, the finite- 
temperature expansion first breaks down as p ^ or 
p — > 1. To more carefully test the validity of the low- 
temperature expansion we also perform the full mean- 
field calculation by numerically solving the nonlinear self- 
consistent equations. The results are displayed along 
with the lowest order calculations in Fig. [5l In order to 
assess the equilibrium stability of the bicontinuous gy- 
roid in the finite temperature phase diagram, we plot the 
grand potential for hexagonal and bcc phases, which are 



the stable phases in the strong segregation regime. The 
full calculation results for lamellar, gyroid, hexagonal, 
and bcc phases are presented in Fig. [B] for ~ 1/10000, 
1/5000, 1/2000 and |V;| = 0.2. Despite the trend sug- 
gested by the lowest-order finite temperature corrections, 
we find that the gyroid structure is not stable at finite 
temperature and the topology of the mean-field phase 
diagram is not qualitatively changed in the low temper- 
ature regime due to the presence of finite-temperature 
corrections. 



V. WEAK-SEGREGATION LIMIT 

Intrigued by near stability of the gyroid phase in the 
finite-temperature expansion, we pursue the mean-field 
clustering phase behavior near to the order-disorder tran- 
sition, where the density is nearly constant with only 
small amplitude modulations. Here, we follow Licbler's 
approach to thermodynamics of block copolymer melts 
[25| near the critical point to derive a Landau theory 
expansion for thermodynamic potential in terms of the 
lattice mean-field theory of clustering. As argued by 
Alexander and McTague, the form of the expansion of 
the free energy in terms of periodic density modulations 
near the liquid-to-solid phase transition takes on a rather 
generic form (ssj . Our task is to derive expressions for 
the coefficients of this expansion in terms of the thermo- 
dynamic parameters of the cluster model: f3, Vq, and 
p. This is most concisely performed by first using the 
density functional theory for the Helmholtz free energy 
(fixed density), which is related to the mean-field theory 
of the action in equation [Tj This yields 

f3F^ J £xS{^) + ^ j d?x' p{^)Vss{^-^)p{^), (23) 

where 5(x) is the lattice-gas expression for the local con- 
tribution to the entropy, eq. ([20|) . We write p(x) = 
po4-(5p(x), where / d?x ^p(x) — and po is constant. To 
Legendre transform to the grand canonical ensemble, we 
add a term — Vppo a-nd minimize with respect to po for 
a fixed density variation. Expanding to quadratic order 
in 5p, we find the grand potential 



(Px 
~V 

d^x 



5p\^) (24) 



V(x) 



The coefficients are defined in terms of the function 
A{p) = pln{p/pm) - [pm - /3)ln(l - p/pm). The mean 
value of density, po, is determined by the solution to the 
equation 



VoPo + A^'Hpo)^Pp, 



(25) 
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where ^^"^ = (PA/dp"'. This allows us to write the co- 
efficients as 

r2(q) = A(2)(p-„) + y(q) 

7 = A(3)(po) 

A = A^'Hpo) (26) 

2{Vo+A(^){po) 

Note that the final, non-local quartic term in the Landau 
expansion of grand potential does not appear in origi- 
nal analysis of Alexander and McTague, nor within the 
mean- field treatment of diblock copolymer melts [2^1 • 

At this point we now consider the grand potential for 
various periodic structures. As before, the function of 
y(q) is to select structures with a periodic modulation 
of density at |q| = g*, to minimize the values of the 
quadratic coupling of the density, r2(q)|(5/5(q)p. Thus, 
we consider density variations which, like iplpc), are plane 
wave superpositions of a given symmetry group. Follow- 
ing a similar analysis for copolymer melts [27j , for a given 
symmetry the density modulation can be expanded in 
terms of the basis functions 



6p{x) =^ (X) 



(27) 



where is an amplitude (not related to Icvcl-sct pa- 
rameter) and 0i(x) is given by a superposition of plane 
waves 



(28) 



. Here, are a family of reciprocal basis vectors corre- 
sponding to the symmetry of a candidate structure and 
the CGi are phase factors satisfying IcgJ = 1- We again 
consider only modes |Gi| = q* . 

It was the original insight of Alexander and McTague, 
that the mean-field behavior of the liquid-to-solid tran- 
sition is dominated by structures that contain a large 
number of triplets for which -|- Gj + Gfc = 0. This 
condition allows the cubic term in cq. ((24l) to be max- 
imally negative, lowering the free energy of certain pe- 
riodic structures.. It is this geometrical feature of the 
reciprocal lattice of the bcc structure (8 triplets) which 
generically favors this structure in the mean-field the- 
ory of the liquid-to-solid transition [sl] as well as in 
the phase diagram of weakly segregated block copolymer 
mehs H. 

More recently it was discovered by Morse and cowork- 
ers [1^ [13 that there was a non-cubic bicontinuous struc- 
ture, the orthorhombic Fddd structure, that was both 
unimodal and contained an abundance of fundamental 
mode triplets that add to zero. Structurally, this Fddd 
structure is related to an orthorhombic distortion of the 
cubic double-gyroid structure with Ia3d symmetry which 
enforces the condition that all 14 wavcvectors to have the 




(a) 




FIG. 7: The mean-field phase diagrams of isotropic, clustering 
particles in the limit of weak-segregation. The phase diagrams 
for two soft-shoulder potentials are shown: — 0.2Vb in (a) 
and |t4| = 0.8Vo in (b). All phase transitions are first-order 
except at the critical point (dark circle) where the system 
may transition continuously from the disordered state to the 
lamellar phase. 



same length. The Fddd phase is constructed from three 
families of modes: Gi = 2tt[±V3x/2± V3y/A± z/4]/ X^, 
G2 = 2n[±V3y/2 ± z/2]/A* and G3 = 27r[±z]/A*. 
Therefore, for this structure we also have three indepen- 
dent amplitudes '^i to consider. 



Using the above Landau expansion, we derive the fol- 
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lowing expressions for the grand potential 



^hex 

f^bcc 
^Fddd 



r2(g*)i'^ + (A/4-,5) -^^ 



r2(g*)*' 



3^3 
4 

3\/6 



7*3 



— A-<5 1 *4 
12 



(29) 



5 *^ 



Vl6 ^ 8 ^ 



_^f2^2^'3 



The mean-field phase diagram near the critical point 
is computed by minimizing over for each candidate 
structure and comparing values to find the structure 
which minimizes the potential. 

For two values of Kk/Vq we show the computed phase 
diagrams in Fig. [T) All phase transitions are predicted 
to be first-order except at the critical point /JK = 4 and 
fJ- = Vq/2 where the system is predicted to pass from 
the uniform, disordered state to the lamellar phase via a 
continuous, second-order phase transition. Interestingly, 
near the onset of order the the same cluster lattices are 
stable: lamellar, hexagonal columnar and bcc lattice of 
spherical clusters. More surprising is the fact the mean- 
field theory predicts a narrow window of a stable bicon- 
tinuous morphology, between the lamellar and columnar 
morphologies. 

Thus, we predict a mean-field phase diagram for clus- 
ter particles that is on the whole rather similar to that 
was originally predicted by Leibler [2^ -and updated by 
Ranjan and Morse [23|-for diblock copolymer melts. In 
that system, beyond the weak-segregation theory, it be- 
comes crucial to consider periodic structures with density 
modulations for which |q| ^ q* . In particular, this is a 
critical ingredient for establishing that the doublc-gyroid 
structure, rather than the orthorhombic Fddd morphol- 
ogy, is by far the most prominent equilibrium bicontin- 
uous structure in copolymer melts [33, 11^. Indeed, it is 
quite reasonable to expect that a more realistic model of 
a soft-shoulder potential that includes values of V{q) for 
"off-peak" wavevectors will predict the stable bicontinu- 
ous structure to be the cubic double-gyroid morphology. 
However, this complication is beyond the scope of our 
present analysis. This shortcoming notwithstanding, our 
present analysis is sufficient to conclude-at least within 
some very narrow region of phase space in the neighbor- 
hood of the critical point-that there is a stable bicontin- 
uous morphology within the mean-field phase diagram of 
isotropic, clustering particles. 

Finally, it is important to note that our predictions 
have focused strictly within the limit of mean-field. It 



is well-known that near to the order-disorder transition 
fiuctuations of the order parameter-long-wavelength ex- 
citations of the ordered structures-as well as fluctuations 
within the disordered play an important role in deter- 
mining phase behavior [3J] . Within the context of block 
copolymer melts, for example, these fluctuations shift 
the thermodynamic onset of periodic order to slightly 
lower temperatures and generally promote the stability 
of the lamellar morphology with its abundance of soft- 
modes nil. It would be straightforward to adapt 
those theoretical methods to the analyses of clustering 
mesophases, and expect the effect of thermal fluctuations 
to play a similar role in the present context. 



VI. DISCUSSION 

We have presented calculations of the three- 
dimensional phase diagram of a lattice gas of particles 
interacting through a soft corona. When the interac- 
tion potential displays a Fourier transform with a sharp, 
negative peak, we find a universal phase diagram in 
both the strong-segregation (low temperature) and weak- 
segregation (near critical) limits. These phase diagrams 
mirror those found in many soft systems, most promi- 
nently for diblock copolymer melts. In the limit of weak- 
segregation, this correspondence becomes more concrete 
as the mean-field model describing the onset of peri- 
odic order has the same generic form for both systems. 
Though our lattice model is unable to distinguish be- 
tween the different core packings leading to these com- 
plex ground states, our results corroborate generally nu- 
merical calculations by Pauschenwein and Kahl [l3, [HI 
of the ground state of particles with a particular soft- 
shoulder potential, the step- function coronas described in 
Eq. ^ . A notable divergence of these analyses occurs at 
very high particle densities, or large chemical potential. 
Our results demonstrate a generic equivalence-even at 
low temperature-between two- and three-dimensionally 
modulated cluster morphologies and their inverse struc- 
tures [see Fig. ([3])]. These inverse morphologies are ab- 
sent from the high-density predictions of refs. [l3l and 
[ll| . but appear in finite-temperature, Monte Carlo simu- 
lations of step-shoulder potentials in two-dimensions [ll] . 

Moreover, we are able to obtain systematic corrections 
to the low temperature behavior and find that these act 
as an effective negative surface tension - thermal fluctu- 
ations prefer to create interface. This occurs because the 
entropy is concentrated along the cluster interfaces. This 
surface tension decreases the gyroid free energy but docs 
not, apparently, stabilize it at temperatures in which our 
low temperature expansion is valid. The low energy of 
the gyroid compared to other bicontinuous phases can 
be attributed to the fact that the trigonometric approxi- 
mation to the gyroid has a smaller mean curvature, and 
is therefore an extremum of the area, compared to other 
symmetry structures. This scenario is similar to a "pack- 
ing frustration" which occurs in block copolymer mor- 
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phologies [3^, [4^. In that case, the propensity for an 
extremal interface between two immiscible blocks is frus- 
trated by the stretching of the polymers. Apparently, 
the double gyroid is the extremal surface in which the 
polymers are the most uniform. In the case of cluster- 
ing particles, this three-dimensional uniformity is related 
to the presence of a single characteristic length scale in 
the interaction potential and the entropy frustrates this 
energetic contribution. 

For clustering interactions in which a single length 
scale dominates, we can also turn these results around to 
give a geometrical interpretation of cluster morphologies 
in the strong-segregation regime. On the one hand, ener- 
getics (roughly) minimize the number of corona overlaps 
while, on the other hand, entropy maximizes interfacial 
area. 

In the weak-segregation regime, coarse-grained ap- 
proaches to colloidal assembly with competing interac- 
tions of various levels of sophistication have predicted 
universal phase diagrams similar to ours [H, 0, [s^ • We 
are aware of one prediction in this regime of a possi- 
ble bicontinuous phase with cubic symmetry [33 |. Our 
results give this universality a compelling geometrical in- 
terpretation at low temperatures in terms of minimizing 
the number of overlaps while simultaneously maximiz- 
ing interfacial area under appropriate constraints. The 
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